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SECTION  1 


INTRODUCTION 


Under  normal  circumstances,  the  density  of  the  earth's  atmosphere 
decreases  monotonically  with  increasing  altitude. ^  This  vertical 
density  gradient,  with  the  addition  of  a  highly  variable  contribution 
due  to  the  water  vapor  component,  is  responsible  for  a  gradient  in  the 
index  of  refraction  at  radio  and  microwave  frequencies.  Typically, 
electromagnetic  waves  in  the  atmosphere  will  bend  downwards,  i.e., 
towards  the  denser  medium.  For  a  generalized  atmosphere  whose  refrac- 
tivity  is  a  function  of  three  dimensions,  elaborate  mathematical  pro¬ 
cedures  are  necessary  to  predict  the  propagation  path  of  radio  waves 
with  reasonable  confidence.  If  the  refractivity  gradient  is  constant 
and  vertical  everywhere,  which  corresponds  to  a  stratified  atmosphere 
in  which  the  refractivity  decreases  linearly  with  altitude,  the  ele¬ 
gant  and  widely  used  effective  earth  radius  method2  (EERM)  may  be 
employed  to  obtain  geometrical  quantities  which  characterize  the  prop¬ 
agation  path,  e.g.,  earth  grazing  angle.  Since  the  EERM  is  predicated 
upon  average  atmospheric  conditions,  it  remains  unsuitable  for  highly 
accurate  prediction  of  propagation  paths  if  detailed  atmospheric  data 
can  be  used  in  a  ray  trace  analysis. 

The  EERM  allows  the  propagation  of  radio  waves  to  be  represented 
by  straight  lines  if  an  appropriate  factor  is  used  to  scale  the 
earth's  radius  as  measured  to  the  terrain  above  which  propagation 
occurs.  When  the  correct  scale  factor  is  used,  slant  range,  grazing 
angle  and  ground  range  are  identical  to  the  true  curved  path  values, 
but  can,  however,  be  reckoned  easily  using  straight  lines.  Although 
the  EERM  is  simple  to  use,  its  accuracy  hinges  upon  the  scale  factor 
used,  and  it  is  limited  to  stratified  atmospheres  which  do  not  exhibit 
excessive  anomalies,  e.g.,  thermal  inversions,  inhomogeneities,  etc. 
Despite  these  reservations,  radar  performance  analysts  routinely  adopt 
a  scale  factor  which  can  vary  from  ~1.1  to  ~1.8,  depending  upon  sea¬ 
son  and  geographic  location.*  The  often  quoted  factor  of  4/3  is  some¬ 
times  used  as  a  very  general  "rule  of  thumb"  by  those  who  lack 
detailed  data  or  wish  to  consider  only  the  general  propagation  envi¬ 
ronment. 


*This 'range*  of  "scale  factors  should  not  be  confused  with  scale  factor 
altitude  dependence  at  a  given  location. 


For  approximate  work,  it  would  be  desirable  to  maintain  the  sim¬ 
plicity  of  the  EERM,  while  extending  the  range  of  the  method  to  higher 
altitudes  for  which  the  refractivity  profile  is  exponential  rather 
than  linear.  The  purpose  of  the  present  work  is  to  show  that  suitable 
scale  factors  can  be  determined  and  employed  to  yield  reasonable 
approximations  for  path  parameters  at  higher  altitudes.  In  general, 
the  scale  factor  is  strongly  dependent  upon  the  source  (transmitter) 
altitude  and  surface  refractivity.  The  above  findings  are  based  on 
numerical  ray  trace  computations  performed  in  the  context  of  a 
semi -empirical  refractivity  profile,  as  described  in  detail  in  the 
next  section. 

It  should  be  emphasized  that  the  EERM  is  no  substitute  for  the 
detailed  atmospheric  modeling  and  ray  trace  techniques  required  for 
very  accurate  propagation  path  prediction. 

In  section  3  results  are  presented  which  include  scale  factor 
versus  source  altitude  curves  and  tables.  Conclusions  appear  in  sec¬ 
tion  4,  while  a  description  of  the  refractivity  model  is  provided  in 
the  appendix. 


SECTION  2 


DESCRIPTION  OF  RAY  TRACE  AND  EFFECTIVE 
EARTH  RADIUS  METHODS 


Although  a  large  body  of  literature  exists  on  the  subject  of 
radio  and  microwave  propagation  in  the  atmosphere,  a  comprehensive 
review  of  this  material  is  unnecessary  to  achieve  the  limited  goals  of 
the  present  work.  Indeed,  the  ray  trace  computations  outlined  below 
require  only  elementary  geometry,  Snell's  law,  and  a  realistic  atmo¬ 
spheric  refractivity  profile.  The  central  purpose  of  this  study  is  to 
investigate  effective  earth  radius  approximations  for  higher  altitude 
paths  and  determine  suitable  scale  factors.  Other  topics  such  as 
attenuation,  electrical  path  length,  and  divergence  are  not 
addressed. 

Consider  a  source  of  radio  or  microwaves  (e.g.,  radar)  at  an 
altitude  h7  above  mean  sea  level  (MSL)  in  a  stratified  atmosphere 
which  lies  above  smooth  spherical  terrain  with  elevation  hs  (MSL), 
as  depicted  in  figure  1.  Irregular  terrain  is  also  permissible  pro¬ 
vided  that  the  terrain  roughness  is  small  relative  to  the  lower  con¬ 
centric  shell  thickness,  4h.  As  illustrated  in  the  figure  (exagger¬ 
ated  for  clarity),  the  spherically  symmetric  atmosphere  is  divided 
into  a  set  of  7  shells  such  that  the  thickness  ratio  of  adjacent 
shells  is  constant.  The  first  shell  thickness  (4h),  and  number  of 
shells,  M,  (M  =  7  in  figure  1)  are  independent  parameters  selected  by 
the  atmospheric  modeler.  These  parameters,  toqether  with  the  total 
thickness  (*17  -hs),  determine  the  adjacent  shell  thickness  ratio, 

C,  through, 


hy  hg 


4h(C#-l  ) 
C  —  1 


(1 ) 


which  can  be  solved  numerically.  The  geometrical  progression  in  shell 
thickness  provides  more  shells  at  the  lower  altitudes,  for  which  the 
refractivity  gradient  is  generally  larger,  to  improve  the  accuracy  of 
the  ray  trace  algorithm.  For  clarity,  only  7  shells  are  illustrated 
in  figure  1.  In  actual  calculations,  as  many  as  100  are  used.  Within 
each  shell  or  stratum  the  refractivity  is  taken  to  be  constant  with  a 
value  specified  by  the  semi -empirical  refractivity  function,  described 
In  the  appendix,  evaluated  at  the  average  shell  altitude.  Therefore, 
in  the  model,  rectilinear  propagation  occurs  within  shells,  while  at 


7 


the  boundaries  Snell's  law  of  refraction  is  invoked  to  account  for  rav 
bendinq.  Ray  tracinq  is  performed  as  follows: 

1.  The  initial  ray  elevation  angle  0M,  source  height  ( h^ ) , 
surface  height  (hs),  number  of  shells  (M)  and  first  shell 
thickness  (Ah)  are  specified. 

2.  The  adiacent  shell  thickness  ratio  (C)  is  obtained  by  numeri¬ 
cal  solution  of  equation  1.  The  thickness  of  the  "mth"  shell 
is  then: 

4hm  =  4hCm_)  (2) 


3.  The  intersection  of  the  ray  with  the  first  shell  boundary, 
i.e.,  the  boundary  between  shell  M  and  M  -1,  is  found  by 
geometry.  The  quantities  of  interest  are  rM,aM  and  #M. 

The  quantity  J~M  -jS  the  length  of  the  segment  of  ray  within 
shell  M  and  is  given  by, 


in  which  Re  is  the  nominal  earth  radius  of  6,373  km.  The  angle 
is  subtended  at  the  earth's  center  by  the  endpoints  of  rM  and  is 
given  by, 


a 


M 


-1 

sin 


rMCOS<^M 

Re+hM-4hCM_1 


(4) 


while  the  incidence  angle  at  the  lower  shell  boundary,  0^,  is  given 

by, 


*M=-aM 


(5) 
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which  Is  a  positive  quantity.  This  is  true  since  only  negative  values 
of  <t> m,  the  initial  elevation  angle,  are  used,  is  always  posi¬ 
tive,  and. 


\K\  >  |“m 


(6) 


from  geometry. 

4.  At  the  first  shell  boundary  the  ray  undergoes  refraction  in 
accordance  with  Snell's  law  (the  interface  is  assumed  to  be 
locally  flat).  Thus,  the  elevation  of  the  ray  at  the  top  of 
the  lower  shell  is  obtained  from. 


in  which  n^  is  the  refractive  index  of  the  upper  shell  and 
njv|_i  is  the  refractive  index  of  the  lower  shell.  Note 
that0M_i,  like  0m,  is  a  negative  (elevation)  angle. 

The  refractive  index  within  a  shell  is  obtained  from  the 
refractlvity  profile  N  [h],  evaluated  at  the  average  shell 
altitude,  for  example; 


nM  —  1  +  N[h] 


h  =  hM 


4hCM  1 
2 


(8) 


5.  Steps  (1)  through  (4)  are  repeated  for  shell  M-l  to  determine 
ctM-1  and  #M-1  >  and  procedure  is  carried 
out  for  the  remaining  shells,  in  sequence. 

In  the  course  of  the  computations  outlined  above,  the  intrashell 
segment  lengths  and  their  corresponding  earth-center  subtended  angles 
are  accummulated.  For  a  large  (~100)  number  of  shells,  the  spatial 
slant  range  from  source  to  surface  is,  to  a  very  good  approximation: 


The  ground  range  from  a  point  on  the  surface  directly  below  the  source 
to  the  terminus  is  given  by. 


M 

Rg  3  (Re+hs)^  am  5  («m  in  radians)  <10) 

m=l 


The  earth  grazing  angle,  6„,  is  simply  the  incidence  angle  at  the 
lower  boundary  of  the  lowest  shell,  l.e.,  shell  no.  1: 

og  =  e]  on 


Thus,  for  given  values  of  source  altitude  ( h^) ,  terrain  alti¬ 
tude  (hs),  and  0m,  the  ray  trace  propagation  algorithm  yields 
Re,  Rg  and  On.  On  the  other  hand,  it  should  be  noted  that  for 
elevation  angles  larger  than  some  critical  value,  which  depends  upon 
h^  and  hs,  some  ray  segment  falls  to  Intersect  a  lower  shell,  and 
rm  as  given  by  equation  (2)  is  complex.  In  this  case,  the  ray 
returns  to  upper  shells  and  diverges  from  the  earth's  surface.  There¬ 
fore,  no  earth  grazing  result  exists  for  elevation  angles  above  a 
critical  value,  which  corresponds  to  the  radio  horizon,  in  the  context 
of  the  shel 1  model . 

The  EERM  is  a  shortcut  but  rigorous  means  of  talcing  into  account 
atmospheric  refraction  at  low  altitudes,  where  the  atmospheric  refrac- 
tivlty  profile  is  (Ideally)  linear.  Under  these  circumstances,  elec¬ 
tromagnetic  waves  can  be  regarded  as  propagating  along  straight  lines 
If  the  earth  radius  Is  properly  scaled,  and  ground  range,  slant  range 
and  grazing  angle  have  correct  values. 

Figure  2  Illustrates,  in  an  exaggera'  manner,  the  EERM  geome¬ 
try-  The  scale  factor,  K ,  is  applied  to  „ne  earth's  radius  as  mea¬ 
sured  to  the  surface  over  which  propagation  occurs.  The  calculation 
of  Rc,  Rg  and  0 g  is  then  analogous  to  that  performed  within  a 
single  shell  for  the  ray  trace  model  (equations  (3),  (4),  and  (5)): 


p 

2 

r 

^|((hs+  Re)  +h  —  hsjsin<p 

—  (h  — hs j 

3 

K(hs^Re)^h-hsJ-(h-hs) 

ia-m.mi 


Re  =  K(Re+  hs)sin  1  - - - - — 

g  L  K(Re+  hs)  J 


(13) 


04) 


For  paths  within  one  kilometer  of  the  surface,  the  correct  scale 
factor  is  a  function  of  the  refractivity  at  the  surface  (Ns),  which 
establishes  the  gradient  (see  appendix],  and  under  ideal  conditions, 
is  given  by  the  convenient  expression:3 


K  =  [l  -  0.04665e0,005577Ns]  1  (15) 


The  often  quoted  K  -  4/3  corresponds  to  Nr  s  301.  Surface 
refractivity  typically  lies  between  200  and  400  N-units. 

Calculation  of  the  ground  range  and  slant  range  to  the  radio  hor 
Izon  is  simple,  since  at  the  horizon  Oa  =  0  and  «  =  -<f>,  i.e.,  the 
triangle  of  figure  2  is  a  right  triangle.  The  ranges  for  propagation 
to  the  horizon  are. 


Rs,horz.  =[(h-hs)[2KRe  +  (2K-l)hs  +  h]J  (16) 


-if  h-hs  "I 

?,horz.  ~  K(Re  +  hs)sec  i  +  — __j 


while  the  corresponding  elevation  angle  at  the  source  is: 


0  *  —  sec 


1  [  1  +  — - — — - -1 

K(Re+  M 


(17) 


(18) 


The  EERM  greatly  simplifies  atmospheric  propagation  geometry  and 
it  provides  a  good  approximation  for  propagation  path  parameters  if 
the  proper  scale  factor  is  used  as  is  discussed  in  the  next  section. 
As  noted  previously,  K  =  4/3  is  sometimes  employed  as  a  general  rule 
of  thumb  --  often  without  sufficient  justification.  In  section  3, 
results  produced  by  the  ray  trace  model  are  used  to  determine  the 
altitude  dependence  of  K  and  the  accuracy  of  the  EERM  for  higher  alti 
tude  paths. 


SECTION  3 


RESULTS:  A  COMPARISON  OF  RAY  TRACE  AND  EERM  PREDICTIONS 


Ray  trace  computations  based  on  the  concentric  shell  atmospheric 
model  have  been  performed  using  an  HP-41CV  system.  Input  variables 
are: 

1.  Source  Altitude  (MSL) 

2.  Terrain  (Surface)-  Altitude  (MSL) 

3.  Ray  (beam)  Initial  Elevation  (or  Depression)  Angle 

4.  Refractivity  at  the  Surface  (Ns) 

5.  Number  of  Shells  (M) 

6.  Thickness  of  Lowest  Shell  (4h) 

From  80  to  100  shells  between  the  source  and  surface  were  suffi¬ 
cient,  as  judged  by  the  observed  rapid  convergence  of  the  ray  trace 
algorithm.  The  variation  of  shell  thickness  was  chosen  to  be  roughly 
in  inverse  proportion  to  the  range  of  the  refractivity  gradient  from 
the  source  to  surface  so  that  the  thinnest  shells  were  always  near  the 
surface. 

Ray  trace  program  outputs  included: 

1.  Spatial  Slant  Range  (equation  9) 

2.  Ground  Range 

3.  Earth  Grazing  Angle 

In  tables  1,  2,  and  3,  ground  ranges  and  grazing  angles  obtained 
from  ray  trace  calculations  are  displayed  for  source  altitudes  of  15 
kft.,  45  kft.,  and  60  kft,  respectively.  All  of  these  cases  are  based 
on  a  terrain  altitude  of  1  kft.  and  surface  refractivity  of  300 
N-units.  As  the  Initial  ray  depression  is  decreased,  the  ground  range 
Increases  and  grazing  angle  decreases  as  expected.  Included  for  com¬ 
parison  in  these  tables  are  the  corresponding  EERM  results  obtained 
with  scale  factors  of  1.209,  1.116  and  1.089  for  tables  1,  2  and  3, 
respectively.  The  EERM  fit  to  the  ray  trace  results  was  enforced  at  a 
ground  range  value  of  80%  of  the  EERM  predicted  radio  horizon  ground 
range.  For  example,  in  table  1  the  EERM  radio  horizon  lies  at  a 
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ground  range  of  256.38  km,  and  the  EERM  scale  factor  was  chosen  to 
produce  agreement  at ~ 200  km  as  can  be  seen  for  the  fifth  table  entry. 
Since  the  EERM  horizon  distance  is  itself  dependent  upon  the  scale 
factor,  several  manual  iterations  were  required  to  obtain  the  desired 
fit.  Any  attempt  to  fit  the  EERM  to  the  ray  trace  results  involves  a 
compromise  between  the  agreement  obtainable  at  high  and  lower  ranges. 
The  procedure  described  gives  very  satisfactory  agreement  for  ranges 
out  to  85%  of  the  radio  horizon,  and  fair  to  bad  agreement  over  the 
remaining  interval.  A  fit  to  ranges  near  the  horizon  would  produce 
poor  agreement  at  intermediate  ranges.  It  is  interesting  to  note  that 
the  ray  trace  and  EERM  horizon  ground  ranges  are  in  relatively  good 
agreement  as  can  be  gleaned  from  the  lowest  entries  in  tables  1,  2, 
and  3. 

Comparison  of  ray  trace  and  EERM  predictions  is  far  more  striking 
when  displayed  in  the  format  of  figures  3,  4,  and  5,  in  which  spatial 
slant  range  is  plotted  against  initial  depression  angle.  For  the  rel¬ 
atively  shallow  depression  angles  spanned  by  the  figures,  little  accu¬ 
racy  would  be  lost  if  "ground  range"  were  substituted  for  "slant 
range"  along  the  ordinate,  since  the  rays  are  roughly  parallel  to  the 
ground.  The  fitted  values  (K  ■  1.209,  figure  3;  1.116,  figure  4; 
1.089,  figure  5)  conform  impressively  to  the  ray  trace  results,  while 
the  "overused"  4/3  value  is  clearly  inferior  in  each  case.  The  K  =  1 
plots  correspond  to  a  line  of  sight  {in  vacuo)  calculation.  Note  that 
K  *  4/3  Is  more  accurate  than  K  =  1  for  the  low  altitude  source 
(figure  3)  while  K  =  1  is  more  accurate  than  K  =  4/3  for  the  high 
altitude  source  (figure  5)  as  expected. 

In  figure  6,  earth  grazing  angle  is  plotted  against  initial 
depression  angle  for  the  45  kft.  source  altitude  case.  Once  again, 
the  superiority  of  the  accuracy  obtained  with  the  fitted  scale  factor 
Is  Immediately  evident. 

These  encouraging  findings  prompted  a  more  thorough  investigation 
of  the  scale  factor's  dependence  on  both  source  altitude  and  surface 
refractivity .  Table  4  contains  ray  trace  fitted  K  values  at  nine 
altitudes  for  each  of  three  N  values:  N  =  200,  300,  and  400 
N-units. 

The  low  altitude  case  (4  kft.)  of  table  4  is  of  particular  inter¬ 
est,  since  the  propagation  path  lies  entirely  within  the  linear  region 
of  the  refractivity  profile,  which,  in  the  semi -empirical  model, 
extends  to  1  km  (3.3  kft.)  above  the  surface.  Thus,  for  this  case  the 
ray  trace  model  fitted  scale  factors  should  corroborate  the  prediction 
of  equation  (15),  which  is  based  on  analytical  work.  Indeed,  the 
tabular  values  1.165  (Ns  =  200),  1.327  (Ns  =  300)  and  1.756  (Ns  =  400) 
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Effective  Earth  Radius  Factor  Versus  Source  Altitude 
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are  In  excellent  agreement  with  the  equation  (15)  derived  values  of 
1.166,  1.331  and  1.767,  respectively  -  a  result  which  underscores  the 
accuracy  of  the  ray  trace  model . 

Figure  7  displays  smooth  curves  of  scale  factor  versus  source 
altitude  for  targets  or  receivers  on  1  kft.  terrain.  These  functions 
were  obtained  via  a  manual  fit  to  the  plotted  reference  points  of 
table  4.  For  Ns  *  400  and  300,  the  scale  factor  decreases  monotoni- 
cally  with  Increasing  source  altitude  as  anticipated,  while  for 
Ne  =  200,  a  relatively  flat  region  appears  between  30  kft.  and  50 
kft.  altitude.  The  unusual  variation  of  K  for  the  latter  case  can  be 
understood  as  a  consequence  of  the  refractivity  model  employed  in  the 
calculations.  The  bending  power  of  the  atmosphere  is  proportional  to 
the  magnitude  of  the  refractivity  gradient,  not  the  absolute  refrac¬ 
tivity.  As  seen  in  figure  A-3,  the  gradient  for  the  Ns  =  200  case 
displays  a  sharp  discontinuity  at  9  km,  which  indicates  that  the  bend¬ 
ing  power  of  the  atmosphere  is  larger  over  a  wide  range  of  higher 
altitudes  than  it  Is  at  lower  altitudes.  This  accounts  for  the  flat 
spot  in  the  Ns  =  200  curve  of  figure  7. 

The  scale  factors  tend  to  converge  for  high  altitude  sources 
since  Ns  Influences  the  refractivity  profile  only  up  to  altitudes  of 
9  km  (see  appendix). 

As  noted  earlier,  the  EERM-ray  trace  agreement  degrades  near  the 
radio  horizon.  This  trend  is  illustrated  in  figures  8  and  9,  which 
display  EERM  ground  range  plotted  against  ray  trace  model  ground 
range  for  common  values  of  depression  angle.  Deviations  from  a  45° 
line  through  the  origin  indicate  discrepancies,  and  it  is  seen  that  in 
the  horizon  region  the  EERM  consistently  predicts  ground  ranges  which 
are  somewhat  longer  than  those  predicted  by  the  ray  trace  model.  This 
result  is  a  consequence  of  the  exponential  refractivity  altitude 
dependence  used  in  the  ray  trace  work.  Bending  in  the  lower  atmo¬ 
sphere  is  accentuated,  causing  ground  range  to  fall  short  of  the  EERM 
predictions.  Note  that  for  the  4  km  source  altitude  case  of  figure  9, 
no  discrepancy  Is  evident,  since  the  model  refractivity  gradient  is 
constant  from  the  surface  to  1  km  (3.3.  kft.). 

The  departure  of  the  EERM  approximation  from  the  ray  trace  pre¬ 
diction  near  the  radio  horizon  for  the  fitted  scale  factors  is  not 
especially  disturbing,  since  no  model  can  be  expected  to  produce  con¬ 
sistently  accurate  predictions  in  this  region,  under  ordinary  circum¬ 
stances  anyway.  In  the  horizon  region  irregularities  in  the  refrac¬ 
tivity  structure  are  magnified  as  a  result  of  low  grazing  angles 
(e.g.,  consider  highway  mirages  seen  at  low  grazing  angles  over  hot 
pavement).  As  a  practical  rule,  the  EERM  predictions  may  be  regarded 
as  good  approximations,  when  atmospheric  irregularities  are  not  large. 


25 


for  ranges  up  to  85%  of  the  horizon  distance  at  higher  altitudes 
(90  kft.  >  h  >  60  kft.)  and  for  ranges  up  to  90%  of  the  horizon 
distance  Tor  Tower  altitudes  (4  kft.  _<  h  <  60  kft.). 

The  results  presented  thus  far  apply  for  a  surface  altitude  of  1 
kft.,  which  was  chosen  as  a  nominal  value.  Table  5  displays  a  set  of 
scale  factors  which  were  determined  for  the  same  conditions  as  those 
of  table  4,  save  for  a  specified  surface  elevation  of  2  kft.  (0.61 
km).  For  high  altitude  sources,  the  scale  factor  changes  very  little 
(compare,  table  4)  while  for  intermediate  and  low-intermediate  source 
altitudes,  changes  are  more  significant  but  still  unremarkable.  For 
the  lowest  altitude  case  (4  km),  no  change  in  scale  factor  is  seen, 
since  in  the  linear  profile  region  the  scale  factor  does  not  depend  on 
the  source-surface  altitude  difference.  In  the  context  of  the  present 
atmospheric  model,  the  change  in  surface  altitude  is  approximately 
equivalent  to  an  equal  and  opposite  change  in  source  altitude.  From 
the  slopes  of  the  curves  of  figure  7,  it  is  then  evident  that  the 
scale  factor  will  vary  more  strongly  for  the  lower  altitude  points 
(save  for  the  lowest)  when  the  surface  altitude  is  changed. 


Table  5 


Effective  Earth  Radius  Factor  Versus  Source  Altitude 
and  Surface  Refractivity  For  2  kft.  Terrain  Elevation 


SECTION  4 


CONCLUSIONS 


The  effective  earth  radius  method  (EERM),  which  is  rigorous  for 
propagation  in  atmospheres  with  a  linear  refractivity  profile  (con¬ 
stant  gradient),  has  been  extended  to  higher  propagation  paths  for 
which  the  refractivity  decreases  exponentially  with  altitude.  The 
extension  is  based  on  ray  trace  computations  which  treat  the  earth's 
atmosphere  as  a  collection  of  concentric  shells  whose  refractive 
indices  are  assigned  in  accordance  with  a  semi -empirical  model.4 

The  popular  4/3  EERM  scale  factor  is  only  appropriate  for  prop¬ 
agation  paths  within  1  km  (3.3  kft.)  of  the  surface  when  the  ground 
level  refractivity  is  near  the  "average"  value  of  300  N-units.  For 
higher  altitude  paths,  significantly  smaller  scale  factors  than  those 
obtained  from  calculations  in  the  low  altitude  region  are  indicated  by 
the  ray  trace  results.  Appropriate  scale  factors  have  been  computed 
for  source  altitudes  from  4  kft.  to  90  kft.,  for  surface  refractivit- 
ies  of  200,  300,  and  400  N-units  (tables  4  and  5).  Interpolation  may 
be  used  to  estimate  scale  factors  for  intermediate  cases. 

As  a  cautionary  note,  it  must  be  emphasized  that  the  use  of  aver¬ 
age  or  idealized  atmospheric  refractivity  models,  while  permitting 
computational  simplicity,  does  not  do  justice  to  the  real  complexity 
and  variability  of  the  atmosphere  -  particularly  the  troposphere. 

Model  calculations  based  on  such  idealized  conditions  should  be  used 
only  when  more  detailed  data  are  unavailable,  or  for  overall  conve¬ 
nience  when  accuracy  is  not  crucial. 

Ground  range  estimates  obtained  from  the  extended  EERM  can  not  be 
expected  to  have  high  accuracy  for  basically  2  reasons: 

1.  The  EERM  is  only  an  approximation  to  the  ray  trace  results 
for  altitudes  above  1  km. 

2.  The  ray  trace  results  themselves  are  based  on  a  model  atmo¬ 
sphere  which  only  reflects  average  conditions. 

If  high  accuracy  is  required,  ray  trace  analysis  must  be  performed  in 
the  context  of  a  realistic  model  atmosphere  tailored  to  the  area  of 
interest. 
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APPENDIX 

SEMI-EMPIRICAL  ATMOSPHERIC  REFRACTIVITY  MODEL 


The  refractivity  of  the  earth's  atmosphere  can,  on  the  average, 
be  represented  by  the  model  described  below,  which  is  taken  directly 
from  the  Radar  Handbook. * 

The  atmosphere  is  assumed  to  be  stratified,  i.e.,  it  is  uniform 
at  a  given  altitude,  with  respect  to  latitude  and  longitude  over  the 
earth's  surface.  The  model  refractivity  is  a  function  of  only  three 
variables: 

1.  Altitude  above  mean  sea  level;  h 

2.  Terrain  elevation  above  mean  sea  level;  hs 

3.  Refractivity  at  the  surface  (h  =  hs);  Ns 

Three  altitude  regions  are  defined  by  the  refractivity  model. 

1.  The  linear  region 

2.  The  surface  refractivity  and  surface  altitude  dependent 
exponential  region 

3.  The  fixed  exponential  region 

Each  of  these  regions  is  described  quantitatively  below. 

A.  The  Linear  Region 

Within  one  kilometer  (3.3  kft.)  of  the  earth's  surface  the 
refractivity  Is  assumed  to  decrease  linearly  with  altitude.  This  is, 
therefore,  the  region  in  which  the  effective  earth  radius  model  is 
rigorous  and  In  which  the  scale  factor  is,  therefore,  independent  of 
altitude.  The  refractivity  is  given  by. 


N  =  N$  +  &N(h  -  h$)  ;  h$  <  h  <  h$  +  1 (km)  <a-i> 


In  which, 


^N  =  -7.32  exp(0. 005577  N$) 


(A-2) 


is  the  magnitude  of  the  refractivity  gradient,  which  is  constant  in 
this  region  for  a  given  surface  refractivity. 

B.  The  Ns  and  hs  Dependent  Exponential  Region 

Between  the  top  of  the  linear  region  (hs  +  1  km)  and  9  km  alti¬ 
tude,  the  refractivity  profile  is  assumed  to  be  exponential.  The 
function  must  be  continuous  at  the  boundary  with  the  linear  region  and 
must  give  N  *  105  at  9  km.  The  continuity  requirement  forces  the 
exponential  decay  constant  to  be  dependent  upon  both  Ns  and  hs. 

The  refractivity  is  given  by, 


N  =  N1  exp[c(h  -  h$  - 


in  which. 


h$  -  8/Ln  \105 


and. 


l)j  ;  h$  +  1  <  h  <  9(km)  (a-3) 


i  W» i 


N1  -  Ns  +  AN 


(A-5) 


The  magnitude  of  the  refractivity  gradient  is  also  exponential  and  is 
given  by. 


BN 

Bh 


CN  =  CN-j 


(A-6) 


and  is  not,  in  general,  continuous  at  the  upper  and  lower  boundaries 
of  this  region. 


C.  The  Fixed  Exponential  Region 


The  refractivity  profile  is  exponential  in  the  third,  uppermost 
region,  and  is  independent  of  surface  conditions.  The  decay  constant 
is  fixed  and  the  refractivity  is  105  N-units  at  the  lower  boundary: 


N  =  105  exp  0.1424(9  -  h)J  ;  9<h(km)  (a-7) 


The  magnitude  of  the  refractivity  gradient  is, 


dN 
5h  “ 


-14.952  exp  jo.  1424(9  -  h)] 


(A- 8) 


Figures  A-l  and  A-2  contain  plots  of  refractivity  versus  altitude 
for  a  surface  elevation  of  1  kft.  (0.3048  km)  and  surface  refractivi- 
tles  of  200,  300  and  400  N-units.  Figure  A-3  contains  the  correspond¬ 
ing  plots  of  the  refractivity  gradient,  in  which  the  discontinuities 
at  1.30  and  9  km  are  clearly  indicated.  Although  discontinuities  of 
this  type  are  unphysical,  any  model  of  the  atmosphere  must  make  com¬ 
promises  somewhere,  and  in  this  case  the  gradient  has  anomalies  for 
the  sake  of  maintaining  the  simplicity  of  the  profile. 
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